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(57) Abstract: The invention concerns a method for focusing 
acoustic waves useful for obtaining an image of a field to be 
observed in a dissipative heterogeneous medium (2, 3) around 
which acoustic transducers (Tl-Tn, T'l-T'm) forming an 
imaging network and a target network. The method consists 
in following a training step during which pulse responses from 
the medium are measured between each transducer (H) of the 
imaging network (5) and several transducers (Tj) of the target 
network (6); deducing therefrom reference signals to be emitted 
by the transducers of the imaging network to produce a focused 
acoustic pulse in each transducer of the target network, then 
cumulatively, in determining reference signals to be emitted to 
focus an acoustic pulse on predetermined points in the medium. 
Said reference signals are stored and used subsequently to 
generate an acoustic image of the medium. 

(57) Abreg£ : Precede de focalisation d'ondes acoustiques utili- 
sable pour obtenir une image d'un champ a observer dans un mi- 
lieu heterogene dissipatif (2, 3) autour duquel on fixe des trans- 
ducteurs acoustiques (Tl-Tn, T' l-T'm) formant un reseau d'ima- 
gerie et un reseau cible. On suit une etape d'apprentissage au 
cours de laquelle on mesure des reponses impulsionnelles du mi- 
lieu entre chaque transducteur (Ti) du reseau d'imagerie (5) et 
plusieurs transducteurs (Tj) du reseau cible (6). On en deduit des 
signaux de reference a emettre par les transducteurs du reseau 
d'imagerie pour produire une impulsion acoustique focalisee en 
chaque transducteur du reseau cible, puis de proche en proche, 
on determine des signaux de reference a Smettre pour focaliser 
une impulsion acoustique sur des points predetermines dans le 
milieu. Ces signaux de reference sont memorises ct utilises cn- 
suite pour generer une image acoustique du milieu. 
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Procede et dispositif non invasif de focalisation d y ondes 
acoustiques . 



La presente invention est relative aux procedes et 
5 dispositif s non invasif s de focalisation d' ondes 
acoustiques, notaminent ultrasonores . 

Plus particulierement, l 1 invention concerne un 
procede non invasif de focalisation d' ondes acoustiques 
dans un milieu heterogene dissipatif comprenant un milieu 
10 sensiblement homogene (par exemple, le cerveau) entoure au 
moins partiellement par une couche aberratrice dissipative 
(par exemple, le crane) qui genere des aberrations dans la 
propagation des ondes acoustiques, les ondes acoustiques 
etant emises depuis l'exterieur de la couche aberratrice et 
15 focalisees dans le milieu sensiblement homogene. 

Les procedes de ce type qui sont couramment 
utilises ne permettent pas d'obtenir une bonne focalisation 
des ondes acoustiques a 1 f interieur "du milieu, et " lorsque 
ces procedes sont utilises dans des applications 
20 'd'imagerie, ils ne permettent done pas d'obtenir une bonne 
resolution et un bon contraste d' image lorsque les 
aberrations de propagation sont importantes, par exemple 
lorsqu'on realise une echographie du cerveau depuis 
l'exterieur du crane. 
25 La presente invention a notamment pour but de 

pallier cet inconvenient. 

A cet effet, selon 1' invention, un procede de 
focalisation du genre en question est caracterise en ce 
qu'il comporte les etapes suivantes : 
30 (a) une etape initiale de positionnement au cours 

de laquelle on fixe un nombre t superieur a 2 de 
transducteurs acoustiques dans des positions predetermines 
a l'exterieur de la couche aberratrice, ces transducteurs 
etant en contact (direct ou indirect) avec ladite couche 
35 aberratrice et formant au moins : 
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un reseau d' imagerie qui regroupe un nombre n 
compris entre 1 et t desdits transducteurs, 

et un reseau cible qui regroupe un nombre m 
compris entre 1 et t desdits transducteurs (ces deux 
5 reseaux peuvent etre entierement distincts, ou cornporter 
certains transducteurs communs, ou encore cornporter chacun 
la totalite des transducteurs susmentionnes ) , 

(b) une etape d' apprentissage comprenant elle-meme 
les sous etapes suivantes : 

10 (bl) une sous-etape d' apprentissage de la 

fpcalisation du reseau d' imagerie sur le reseau cible, 
sous-etape au cours de laquelle : 

(bll) on determine des reponses 

impulsionnelles hri(t) du milieu heterogene dissipatif, 

15 respectivement entre chaque transducteur i du reseau 
d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible (cette" determination peut 
etre faite par mesure directe si l'on fait emettre des 

20 impulsions acoustiques par les transducteurs du reseau 
cible, ou le cas echeant par mesure et calcul si l'on fait 
emettre des signaux acoustiques autres que des impulsions 
par les transducteurs du reseau cible, les valeurs ainsi 
mesurees et/ou calculees pouvant ensuite etre le cas 

25 echeant corrigees par repropagation numerique pour simuler 
des transducteurs situes directement au contact de la 
couche aberratrice si les transducteurs ne sont pas au 
contact direct de la couche aberratrice) , ces reponses 
impulsionnelles etant memorisees sous forme numerique avec 

30 un certain echantillonnage temporel qui determine un nombre 
p de composantes f requentielles de la reponse 
impulsionnelle, de frequences respectives ok, i etant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d ! imagerie, r etant un indice compris entre 1 et m 

35 qui designe un point de focalisation correspondant a un 
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transducteur du reseau cible et k etant un indice compris 
entre 1 efc p qui designe une composante f requentielle , 

(bl2) a partir de ces reponses 
impulsionnelles, on calcule, pour chaque point de 
5 focalisation r correspondant a un transducteur du reseau 
cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la paroi 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
10 les differents transducteurs i du reseau d' imagerie 
genererait un signal predetermine (par exemple une 
impulsion acoustique) focalise sur le point de focalisation 
r, 

(b2) une sous-etape de focalisation en un nombre 
15 R de points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 
et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation - q, en s' eloignant pas 
a pas des points de focalisation lam correspondant aux 

0 transducteurs du reseau cible, des signaux de reference 
e'i(t,q> a faire emettre par les differents transducteurs i 
du reseau d' imagerie pour generer une impulsion focalisee 
sur ledit point de focalisation q, les signaux de reference 
e'i(t,q) etant determines pour chaque point de focalisation 

5 q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 
pour i allant de 1 a n, est calculee a partir d'au moins un 
signal de reference e'i(t,q0), qO etant l ! indice d'au moins 
un point de focalisation proche du point de focalisation q 

3 pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
ondes acoustiques dans le. milieu sensiblement homogene, 

(b22) on fait emettre par les transducteurs 
du reseau d' imagerie, par iterations, les estimations 

3 precedemment obtenues des signaux de reference " e'i(t,q) , 
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puis on capte avec les memes transducteurs des signaux 
Sj,(t,q) retrodif fuses par le milieu heterogene dissipatif, 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q> de la maniere suivante : 

5 e\{t) a f (f/).e|(/-r ; .(7)) 

ou les valeurs (^(q) et Ti(q) sont un facteur d' amplitude 
et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses, 
lesdites iterations etant arretees lorsque le critere C 

10 atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 

15 q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i(t,q), i 
etant un indice compris entre 1 et n designant un 
transducteur du reseau d ! imagerie. 

Grace a ces dispositions, on s ! affranchit des 

20 aberrations de propagation des ondes acoustiques dans le 
milieu heterogene dissipatif, et on obtient une 
focalisation tres precise qui peut notamment permettre 
d'obtenir une image echographique fidele et precise d'un 
champ a observer a travers la couche aberratrice par 

25 retrodif fusion, lorsqu'on emet successivement des ondes 
acoustiques focalisees sur differents points du champ a 
observer et que l'on capte les ondes acoustiques 
retrodif fusees . 

Cette focalisation precise peut egalement etre 

30 utilisee dans d'autres applications que 1' echographie , 
notamment : 

l'imagerie doppler couleur, 

les methodes d' imagerie par elastographie , 
telles que celle decrite dans le document WO-A-00/55 616, 
35 - les methodes d' imagerie non lineaire ("harmonic 
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imaging" ) , 

- les methodes de traitement par destruction 
localisee d'une partie du milieu heterogene dissipatif, 
notamment par hyperthermie, 

- les methodes de mesure de parametres 
d' absorption optique des tissus avec activation par 
ultrasons, etc. 

Dans des modes de realisation preferes de 
1' invention, on peut eventuellement avoir recours en outre 
a l'une et/ou a 1' autre des dispositions suivantes : 

- au cours de la sous-etape (bll), lorsqu'au moins 
certains transducteurs (du reseau cible et/ou du reseau 
d'imagerie) sont en contact avec un milieu homogene 
intermediaire (par exemple un gel) lui-meme en contact avec 

15 la couche aberratrice, on corrige les reponses 
impulsionnelles hri(t) par repropagation numerique pour 
simuler des transducteurs situes directement au contact de 
la couche aberratrice ; 

- la sous-etape (bl2) comporte elle-meme les sous- 
20 etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(cok) = [Hri (cok) ] , i allant de 1 a n et r allant de 1 a m, 
ou Hri(cok) est la valeur, a la frequence cok, de la 
transformer de Fourier de la reponse impulsionnelle hri(t), 
25 (bl22) on determine pour chaque point de 

focalisation r correspondant a un transducteur du reseau 
cible, n composantes Ei(cok,r), i variant entre 1 et - n, 
telles que F (cok, r) =H (cok) . E (cok, r) , ou E(cok,r> = (Ei(cok,r)] 
est un vecteur a n composantes, F(cok,r) est un vecteur a m 
30 composantes Fl(cok,r), 1 variant entre 1 et m, ces m 
composantes Fl(cok,r) correspondant a une focalisation 
souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant a un transducteur du 
reseau cible, 
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(bl23) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux ternporels e (t, r> = [ei ( t, r) ] , 

/' 

i variant entre 1 et n, ou e,{i i r)=Y,Ei{a)kj'). e tt ' ' en notation 

5 cornplexe, ces signaux ei(t,r) etant adaptes pour que 
1' emission de ces derniers respectivement par les 
differents transducteurs i du reseau d'imagerie genere une 
impulsion acoustique focalisee sur le point de focalisation 
r du reseau cible, 

10 . (bl24) une sous-etape de correction des 

aberrations generees par la couche aberratrice entre le 
milieu sensiblement homogene et chaque transducteur cible 
r, ces aberrations etant estimees sur la base des mesures 
precedemment effectuees, les aberrations ainsi estimees 

15 etant utilisees pour calculer lesdits signaux ternporels de 
reference e'i(t,r) ; 

- au cours de la sous-etape (bl22) on calcule p 
matrices H"" 1 (cok) , respectivement par regularisation et 
inversion des matrices de transfert H(cok) / et pour chaque 

20 transducteur r du reseau cible, on calcule le vecteur 
E(cok,r) par la formule : 

Eto^rHf-TMcok) . F(cok, j) ; 

- au cours de l ! etape (bl22), les composantes 
Fl(cok,r) du vecteur F(o)k,r) correspondant a la repartition 

25 spatiale du champ desire a la frequence cok, sont egales a 
0 pour l^r et egale a 1 pour l=r ; 

- au cours de la sous-etape (bl24), on assimile la 
paroi aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du reseau 

30 cible, a un filtre a reponse impulsionnelle finie, defini a 
chaque frequence cok par une amplitude Gr{a)k) et une phase 
(|> r (c0k), la sous-etape (bl24) comportant elle-meme les sous- 
etapes suivantes : 
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(bl241) on calcule, pour chaque frequence cok, 
1' amplitude Gr(cok) et la phase <|> r (cok>, a partir soit des 
signaux ei(t,r), soit des vecteurs E(wk, r) , 

(bl242) on calcule p matrices de transfert 
5 corrigees H' (cok) = [H' ji (cok) J , ou 



du reseau cible, n composantes E'i(c»k,r), i variant entre 
1 et n, telles que F(cok,r)=H' (cok) .E' (a)k,r), ou E' (cok, r ) - 

10 [Ei(cok,r)J est un vecteur a n composantes, F(cok,r) est un 
vecteur a m composantes Fl(cok,r), 1 variant entre 1 et m, 
ces m composantes Fl(a>k,r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
<ok sur le point de focalisation r correspondant a un 

15 transducteur du reseau cible, 



focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels de reference 
e' (t,r) = [e'i(t,r>], i variant entre 1 et n, ou 
p 

20 e'i(t i r)=^Fi(o)k,r).g ' en notation complexe ; 



Hp (a> k ) = H n {co k ) . — — - e ~ 



(bl243) on determine pour chaque transducteur r 



(bl244) on en dediiit, pour chaque point de 



1' 



- au cours de la sous-etape (bl241) , on . calcule 
amplitude Gr(o)k) et la phase <J> C (cok) comme suit : 




25 




ou 



' Ei* est la valeur complexe conjuguee de Ei, 
et At(rO,r,i)={d(rO,i>-d(r, i> > /c, d(r, i) etant la 
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distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO ; 

la sous-etape (bl2) comporte elle-meme les sous- 
5 etapes suivantes : 

(bl21) on determine p matrices de transfert 
H (cok) = [Hri (cok) ] , i allant de 1 a n et r allant de 1 a m, 
ou Hri (cok) est la valeur, a la frequence cok, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 

10 (bl22' ) on corrige les matrices de transfert 

H'(cok) pour s'affranchir des aberrations generees par la 
paroi aberratrice au voisinage de chaque point de 
focalisation r, cette correction etant effectuee a partir 
des reponses impulsionnelles hri (t) precedemment 

15 determinees, et on obtient ainsi des matrices de transfert 
corrigees H' (cok) , 

(bl23') on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, n composantes E'i(cok,r), i variant entre 1 et n, 

20 telles que F (cok, r ) =H' (cok) .E' (cok,r), ou E' (cok, r ) = 
[E'i(tok,r)] est un vecteur a n composantes, F(cok,r) est un 
vecteur a m composantes Fl(cok,r), 1 variant entre 1 et m, 
ces m composantes Fl (cok, r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 

25 cok sur le point de focalisation r correspondant a un 
transducteur du reseau cible, 

(bl24') on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels 

30 e' (t, r) -[e' i (t, r) ] , i variant entre 1 et n, ou 

cV(/,/*)=5]£ l /(^%/')e'^' en notation complete, les signaux 

e'i(t,r) etant lesdits signaux de reference ; 

au cours de la sous-etape (bl23 7 > on calcule p 
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matrices H'->(a>k), respectivement par regularisation et 
inversion des matrices de transfert H' (o>k) , et pour chaque 
transducteur r du reseau cible, on calcule le vecteur 
E'(cok,r) par la formule : 

E' (cok, r) «H' -1 (oak) . F (cok, j ) ; 

- au cours de 1'etape (5123'), les composantes 
Fl(cok,r) du vecteur F(a>k,r> correspondant a la repartition 
spatiale du champ desire a la frequence »k, sont egales a 
0 pour l*r et egale a 1 pour l=r ; 

- au cours de la sous-etape (bl22'), on assimile 
la paroi aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, a un filtre a reponse impulsionnelle finie, defini a 
chaque frequence cok par une amplitude Gr(cok) et une phase 
♦r«Dk), la sous-etape (bl22' ) comportant elle-meme les 
sous-etapes suivantes : 

(bl22'l) on calcule, pour chaque frequence 
ok, 1' amplitude Gr(cok) et la phase ^(cuk), a partir des 
reponses impulsionnelles precedemment determines, 

(bl22'2) on calcule p matrices de transfert 
corrigees H' (cok) = [H' ji (cok) } , ou 

au cours de la sous-etape (bl22'l) on calcule, 
pour chaque frequence cok, 1' amplitude Gr(cok) et la phase 
25 (|> r (cok), de la maniere suivante : 



15 



20 



jy j H^.i{m).H,i).t\cok 
V fci 
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H*ri designe la valeur complex e conjuguee de Hri, 
et AT(rO,r,i) = (d(rO,i)-d(r,i) )/c, d(r,i) etant la 
distance entre le transducteur i et . le point de 
focalisation r, et d(rO,i) etant la distance entre le 
5 transducteur i et un point de focalisation particulier rO ; 

- au cours de l'etape (c) , on fait suivre la sous- 
etape (cl) par les sous-etapes suivantes : 

(c2) on fait capter par lesdits t ransducteurs du 
reseau d' imagerie des signaux si(t) retrodif fuses par le 
10 milieu heterogene dissipatif, 

(c3) on convolue le signal de reference emis par 
chaque transducteur du reseau d' imagerie avec le signal 
retrodiffuse capte par ce transducteur, 

(c4) puis on somme les produits de convolution 
15 ainsi obtenus, 

l'etape (c) £tant renouvelee pour une pluralite de points 
situes dans le milieu sensiblement homogene ; 

- au cours de la sous-etape (b21) , la premiere 
estimation de chaque signal de reference est e'i(t,q) = 

20 e' i (ts+0i (q) , qO) pour chaque point de focalisation q, qO 
etant I'indice d'un point de focalisation proche du point q 
pour lequel le signal de reference a deja ete determine, 
0i(q) etant un retard egal a une valeur 5i(q)/c, ou c est 
la celerite moyenne des ondes acoustiques dans le milieu, 

25 et 5i(q) est egal a une difference entre d'une part, une 
distance entre le transducteur i du reseau d' imagerie et le 
point de focalisation qO, et d' autre part, une distance 
entre le transducteur i du reseau d' imagerie et le point de 
focalisation q, 

3 0 - a u cours de la sous-etape (b2>, lorsqu'au moins 

certains transducteurs d 7 indice v du reseau d' imagerie ne 
sont pas directement au contact de la couche aberratrice, 
on corrige les signaux e\-(t,g) correspondants par 
repropagation numerique pour simuler des transducteurs 

35 places en contact direct avec la couche aberratrice ; 
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- au cours de la sous-etape (b22), on recherche 
les valeurs ct x (q) et x,_{q) pour maximiser le critere de 
coherence C suivant : 



C = 







2 


< 




> 









ou 



5 - ft ('.?) = MO®*/ ('.?)# ® representant l'operation de 

convolution, 

et <> represente une moyenne temporelle ; 

- au cours de la sous-etape (b22), les valeurs 
ti(q) sont calculees en maximisant une fonction 

10 d' intercorrelation, pour des transducteurs voisins du 
reseau d'imagerie, des signaux g 1 (t,q) et g i+ i(t,q) ; 

- au cours de la sous-etape (b22>, les valeurs 
cti(q) sont calculees de maniere a egaliser sur l'indice i 

" 1 '"amplitude" maximale des roncf "ions ~q x (t , qj ' / 
.5 - au cours de la sous-etape (b22), les valeurs 

cci(q) et xi(q) sont calculees en realisant une 
intercorrelation, pour des transducteurs voisins du reseau 
d'imagerie, des signaux g x (t,q) et g i+1 (t,q> / 

- au cours de la sous-etape (b22), les valeurs 
0 a x (q) et Ti(q) sont calculees de maniere a egaliser sur 

l'indice i 1' amplitude maximale des fonctions gi(t,q) ; 

- la sous-etape (b22) relative a chaque point de 
focalisation q est realisee immediatement apres la sous- 
etape (b21) relative au rnerne point de focalisation q ; 

5 le milieu heterogene dissipatif est . constitue 

par le cerveau entoure par le crane ; 

le reseau d'imagerie et le reseau cible sont 

deux reseaux distincts disposes de part et d' autre du 

milieu heterogene dissipatif ; 
0. ' tous les transducteurs appartiennent a la fois 

au reseau d'imagerie et au reseau cible ; 
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les ondes acoustiques sont des ondes 
ultrasonores . 

Par ailleurs, 1' invention a egalement pour objet un 
dispositif adapte pour mettre en ceuvre le procede defini 
5 ci-dessus. 

D ! autres caracteristiques et avantages de 
1' invention apparaitront au cours de la description 
suivante d'une de ses formes de realisation, donnee a titre 
d' exemple non limitatif, en regard du dessin joint. 

10 Sur le dessin, la figure unique represente un 

dispositif d'imagerie ultrasonore selon une forme de 
realisation de l 1 invention. 

Le dispositif 1 d'imagerie ultrasonore represente 
sur le dessin est adapte pour realiser une image 

15 echographique par ultrasons du cerveau 2 d'un patient (a 
des frequences par exemple de l'ordre de 1 a 3 MHz), depuis 
l'exterieur du crane 3, le cerveau 2 constituant un milieu 
sensiblement homogene pour la propagation - des ondes 
acoustiques et le crane 3 constituant une couche 

20 aberratrice dissipative, de sorte que 1' ensemble de la 
bolte cranienne 2,3 constitue un milieu heterogene 
dissipatif. 

En variante, 1' invention serait applicable 
notamment : 

25 a 1 1 imagerie acoustique de tout autre milieu 

heterogene dissipatif non homogene comprenant un milieu 
sensiblement homogene entoure par une couche dissipative 
relativement peu epaisse generant des aberrations dans la 
propagation des ondes ultrasonores, 

30 ou a tout autre procede impliquant au moins une 

focalisation a 1' emission dans un tel milieu. 

Dans 1' exemple represents sur le dessin, le 
dispositif d'imagerie 1 comporte un micro-ordinateur 4, ou 
tout autre dispositif de commande et/ou de visualisation 

35 des images ultrasonores, ce micro-ordinateur comportant 
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classiquement un clavier 4a associe eventuellement a 
d'autres interfaces de commande et un ecran 4b permettant 
de visualiser les images du cerveau 2. 

Par ailleurs, le dispositif d'imagerie 1 comporte 
5 deux reseaux 5, 6 de transducteurs ultrasons Tl, T2... Tn et 
T'l, T'2... T'm formant par exemple deux barrettes lineaires 
de transducteurs qui sont disposees de part et d' autre du 
crane 3 de 1 ' utilisateur , dans des positions geometriques 
predeterminees l'une par rapport a 1' autre, chaque barrette 
10 de transducteur 5, 6 etant mise en contact avec le crane 3 
par 1 ' intermediaire d'une couche 7 de gel ou similaire. 

Les differents transducteurs Tl, T2... Tn et T'l, 
T'2... T'm peuvent etre commandes directement par le micro- 
ordinateur 4, ou de preference par une unite centrale 
15 electronique CPU contenue par exemple dans une baie 
electronique 8 et elle-meme commandee par le micro- 
ordinateur 4 . 

Avantageusement, chacun des transducteurs Tl, T2... - 
Tn, T'l, T'2, T'm est relie a un echantillonneur, 
20 respectivement El, E2... En, E'l, E'2, E'm, et chaque 
echantillonneur est lui-meme relie a une memoire, 
respectivement Ml, M2... Mm, M'l, M'2... M'm et a une unite 
centrale CI, C2, ...Cm, C 1, C 2, ... Cm. Ces memoires et ces 
unites centrales sont a leur tour relies, directement ou 
15 indirect ement, a 1' unite centrale CPU susmentionnee, qui 
est par ailleurs reliee au moins a une memoire centrale M. 

Le dispositif qui vient d'etre decrit fonctionne 
comme suit. 

Initialement, les deux reseaux de transducteurs 5, 
;0 6 sont fixes de part et d' autre du crane 3 du patient, dans 
lesdites positions predeterminees. A cet effet, les reseaux 
de transducteurs 5, 6, dits respectivement reseau 
d'imagerie et reseau cible, peuvent etre portes par un 
support rigide tel qu'un casque (non represents) dispose 
5 autour de la tete du patient. 
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Ensuite, le dispositif suit une etape 
d ' apprentissage de quelques minutes (avantageusement 1 a 3 
minutes) permettant de prendre en compte 1' ensemble des 
aberrations de propagation dues au caractere non homogene 
5 du milieu dissipatif forme par le crane 3 et le cerveau 2. 

Au cours de cette etape d* apprentissage on fait 
d'abord emettre successivement par chacun des transducteurs 
Tl, T2... Ti,... Tn du reseau d'imagerie 5, une impulsion 
acoustique, et pour chaque impulsion emise par 1 1 un des 

10 transducteurs Ti du reseau d'imagerie, on enregistre le 
signal capte par les transducteurs T'l, T'L.. T'r,... T'm du 
reseau cible 6, c'est-a-dire la reponse impulsionnelle 
hri(t) du milieu heterogene dissipatif entre le 
transducteur i considere du reseau d'imagerie 5 et chaque 

15 transducteur j du reseau cible 6. 

Chaque reponse impulsionnelle hri(t) est 
enregistree sous forme numerique avec un certain 
echantillonnage temporel qui determine un certain nombre p 
de composantes f requentielles monochromatiques de la 

20 reponse impulsionnelle, correspondant chacune a une 
frequence a>k, k etant un indice compris entre 1 et p. 

Dans le cas envisage ici, ou au moins certains 
transducteurs du reseau cible et/ou du reseau d'imagerie ne 
sont pas directement au contact avec la couche aberratrice 

25 3, on corrige les reponses impulsionnelles pour simuler des 
transducteurs virtuels disposes au contact de ladite couche 
aberratrice. La position de la couche par rapport aux 
transducteurs peut eventuellement etre obtenue par imagerie 
conventionnelle (echographie ultrasonore, scanner X, IRM, 

30 etc.). Les reponses impulsionnelles corrigees sont 
calculees par un algorithme de repropaga tion numerique 
connu, decrit notamment dans les articles suivants : 

"Ultrasonic beam steering through inhomogeneous 
layers with a time reversal mirror" , C.DORME, M. FINK, IEEE 

35 Transactions Ultrasonics , Ferroelectric and Frequency 
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Control, 43 (1), janvier 1996, p 167-175, 

"Focusing and steering through absorbing and 
aberrating layers : Application to ultrasonic propagation 
through the skull" Journal of Acoustical m Society of 
5 America, 103 (5), Mai 1998, p. 2403-2410, 

et "Propagation and backpropagation for 
ultrasonic wavefront design" Liu, D.-L., and Waag, R. C. 
IEEE Trans, on Ultras. Ferro. and Freq. Contr. 44(1):1-13 
(1997) . 

10 Dans ce qui suit, hri(t) denommera done les reponses 

impulsionnelles pour des elements (reels ou virtuels) 
situes contre la couche aberratrice. De plus, les elements 
virtuels ou reels situes contre la couche aberratrice 3 
seront appeles ci-apres "points de f ocalisation" d'indice r 

15 compris entre 1 et m. 

Lorsque les transducteurs du reseau d'imagerie 5 
emettent des signaux acoustiques e A (t>, ces signaux 
gdnerent au niveau des transducteurs r du rfeseau cible 6 
des signaux acoustiques fr(t) s'exprimant corume suit : 

20 /r(0=Shri(t)®ei(t), 

ou ® represente l'operateur de convolution temporelle. 

Apres transformee de Fourier, cette equation 
devient : 

F(o)k) = H(cok) .E(cok) , ou : 

25 - H(ook) est la matrice de transfert, de taille 

m*n, entre les transducteurs Ti du reseau d'imagerie et les 
transducteurs Tr du reseau cible : les composantes Hri(cak) 
de cette matrice sont les composantes des transformers de 
Fourier des reponses impulsionnelles hri(t) a la frequence 

30 ok, 

* * - E(cok} est un vecteur dont les composantes 

Ei(a)k) sont les composantes de la transformee de Fourier 
des signaux e A (t) susmentionnes a la frequence cuk, 
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et F(a>k) est un vecteur dont les composantes 
F-^cok) sont les composantes de la transformee de Fourrier 
des signaux f^(t) susmentionnes a la frequence (uk. 

Par inversion de chaque matrice de transfert H(cok), 
5 on peut done determiner le vecteur E(a)k,j> qui est adapte 
pour generer au niveau du point de focalisation r 
correspondant au transducteur T'r du reseau cible, un 
vecteur F(cok,j) dont toutes les composantes se rapprochent 
au mieux de l'objectif initialement fixe (de preference 
10 toutes egales a zero, sauf la composante d'indice j 
correspondant au transducteur T'j, qui est egale a 1 
lorsqu'on veut emettre une impulsion acoustique au niveau 
du point de focalisation r) , grace a la relation : 

E((ok,j) = H- x (G>k) .F(03k, j), 
15 oil H^ 1 (cok) est la matrice inverse de H(cok). 

IT 1 (ok) peut etre calculee par exemple par 
decomposition en valeurs singulieres, ceci permettant une 
regularisation de 1' inversion de la matrice H(cok). 

Puis on determine, par transformee de Fourier 
20 inverse des differentes composantes Ei(GL>k,j) du vecteur 
E(cok,j), les differents signaux de reference ei(t,j) qui, 
lorsqu ! ils sont emis par les differents transducteurs Ti du 
reseau d'imagerie 5, sont adaptes pour focaliser une 
impulsion acoustique (ou le cas echeant un autre signal 
25 acoustique) au niveau du point de focalisation r. On a done 
realise une focalisation du reseau d' imagerie 5 sur chaque 
transducteur du reseau cible 6, par filtre inverse spatio- 
temporel . 

L' unite centrale CPU suit ensuite un processus 
30 d' apprentissage des aberrations dues a la paroi du crane 3 
au niveau du reseau cible. 

Au cours de ce processus, ces aberrations sont 
considerees comme un filtre a reponse impulsionnelle finie. 

Dans le domaine de Fourier, ce filtre est defini a 
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. chaque frequence cok une amplitude Gr(cok) et une phase 

<|>r(C0k) - 

Pour calculer ces coefficients, on compare la phase 
et. 1' amplitude de 1' ensemble des vecteurs Er. A cet effet, 
5 on commence par eliminer les dephasages introduits par les 
differences de marche entre les transducteurs d'imagerie Ti 
et les differents points de focalisation indices r. Geci 
revient a choisir un point de focalisation particulier rO, 
et introduire pour les autres un dephasage lineaire avec la 

10 pulsation : exp (- jAt (rO, r, i) o) avec At (rO, r, i) - (d (rO, i) - 
d(r,i))/c ou d(r,i) est la distance entre le transducteur i 
et le point de focalisation r et c est la vitesse moyenne 
des ondes acoustiques dans le milieu a imager, * en 
1' occurrence le cerveau 2. 

15 Cette correction effectuee, les differences 

d' amplitude et de phase entre les vecteurs Er sont 
attributes a la couche aberratrice 3 situees contre le 
reseau cible. On calcule alors pour chaque point de 
focalisation r, le facteur de gain Gj (cok) et le facteur de 
20 phase (o>k) : 

A^E{c»k 7 rQyEtXc)k s rt) 



ou Ei* est la valeur complexe conjuguee de Ei. 

25 Ces couples iGj(ok) , $ 5 (©k) ) correspondent au 

facteur d' attenuation et au dephasage relatifs introduits a 
chaque frequence par la portion de couche aberratrice 3 
situee contre le point de focalisation r. lis caracterisent 
done finalement les aberrations introduites par la portion 

30 de couche aberratrice situee contre le reseau cible. 

On elimine ensuite les aberrations introduites par 
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la couche aberratrice 3 situee contre les transducteurs 
cibles, dans 1' ensemble des p matrices H (cok) = t Hj i (ok) ] , 
def inies ci-dessus . 

Pour cela, on calcule un nouvel ensemble de 
5 matrices de transfert H' (cok) = [H 7 ji (ok) ] caracterisant la 
propagation entre le reseau d'imagerie et le reseau cible 
dans un milieu virtuel pour lequel seules subsistent les 
aberrations situees contre le reseau d'imagerie : 

Gj{6) k ) 

.0 On determine alors pour chaque transducteur r du 

reseau cible, n composantes E'i(cok,r), i variant entre 1 
et n, telles que F (cok, r ) =H' (cok) . E' (cok, r ) , ou E'(cok,r) - 
{E' i (cok, r) ] est un vecteur a n composantes, F(cok,r) est un 
vecteur a m composantes Fl(cok,r), 1 variant entre 1 et m, 

.5 ces m composantes Fl(cok,r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant a un 
transducteur du reseau cible. 

On en deduit, pour chaque point de focalisation r 

0 correspondant a un transducteur du reseau cible, un vecteur 
de n signaux temporels de reference e' (t , r ) = [e' i ( t , r) ] , i 

p 

variant entre 1 et n, ou e'f(ts)^^Fi(cvkj').g MJ en notation 

m 

complexe . 

Ces signaux de reference e'i(t,j) sent adaptes pour 
5 que 1' emission de ces derniers respectivement par les 
differents transducteurs i du reseau d'imagerie genere une 
impulsion acoustique focalisee sur le transducteur j du 
reseau cible en 1' absence de la couche aberratrice situee 
contre le reseau cible. 
0 On notera qu' en variante, les signaux de reference 

pourraient etre determines de la facon suivante apres la 
determination des reponses impulsionnelles hri(t) et des p 
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matrices de transfert H(cok) : 

on corrige les matrices de transfert H(cok) pour 
s'affranchir des aberrations generees par la paroi 
aberratrice 3 au voisinage de chaque point de focalisation 
5 r, . cette correction etant effectuee a partir des reponses 
impulsionnelles hri(t) precedemment determinees, et on 
obtient ainsi des matrices de transfert corrigees H' (cok) , 

on determine par inversion des matrices H' (cok) , 
pour chaque point de focalisation r correspondant a un 
10 transducteur du reseau cible, n composantes E'i(cok,r), i 
variant 'entre 1 et n, telles "que F (cok, r) =H' (cok) . E' (cok, r) , 
ou E' (<ok,r) = [E / i(Q)k / r)J est un vecteur a n composantes, 
F(cok,r) est un vecteur a m composantes Fl(a>k, r), 1 variant 
entre 1 et m, ces m composantes Fl(cok,r) correspondant a 
15 une focalisation souhaitee des ondes acoustiques a la 
frequence cok sur le point de focalisation r correspondant 
a un transducteur du reseau cible > 

et on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
20 cible, un vecteur de n signaux temporels e'(t,r) = 

i> 

[e'i(t,r)), i variant entre 1 et n, ou e?i{ts)='£E l iXa)k,r).g t ' en 

notation complexe,. les signaux e'i(t,r) etant lesdits 
signaux de reference. 

Avantageusement, lors du calcul des matrices 
25 H' (cok) , on assimile la paroi aberratrice, au voisinage de 
chaque point de focalisation r correspondant a . un 
transducteur du reseau cible, a un filtre a reponse 
impulsionnelle finie, defini a chaque frequence cok par une 
amplitude Gr(cok) et une phase <f> r (G)k) calcules comme suit : 



30 CrC^)"— 
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^r(a>i0=-^(arg(^ , ou : 

n i»i 

H*ri designe la valeur complexe conjuguee de Hri, 
et At(r0,r,i)=(d(rO,i)-d(r,i) )/c, d(r,i> etant la 
5 distance entre le transducteur i et le point de 
focalisation r, et d(r0,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

On calcule ensuite p matrices de transfert 
corrigees H' (Q)k) s =[H , ji(cDk>] / ou 

10 H) i (o) k ) = H Ji (co k ).—^ e -J*M) f qui S ervent a determiner les 

Gj(o) k ) 

vecteurs E'i(cok) comme explicite precedemment et done les 
differents signaux de reference e'i(t,r), r allant de 1 a 
m. 

L' unite centrale CPU procede ensuite a un 

15 apprentissage de focalisation en un nombre R de points de 
focalisation predetermines situes dans le cerveau 2, 
d 1 indices q compris entre m+1 et m+R, cette sous-etape 
consistant a determiner pour chacun de ces points de 
focalisation q, en s'eloignant pas a pas des transducteurs 

20 du reseau cible, des signaux de reference e'i(t,q) a faire 
emettre par les differents transducteurs du reseau 
d' imagerie pour generer une impulsion focalisee sur ledit 
point de focalisation q. 

Les signaux de reference e'i(t,q) sont initialement 

25 determines, pour chaque nouveau point de focalisation q, 
sous la forme e'i(t,q) = e x i (t+0i (q) , qO ) pour chaque point 
de focalisation q, qO etant 1 ' indice d ' un point de 
focalisation proche du point q pour lequel le signal de 
reference a deja ete determine, le retard 8i(q) etant 

30 initialement egal a une valeur 5i(q)/c, ou c est la valeur 
moyenne de celerite des ondes acoustiques dans le milieu, 
et 5i(q) est egal a une difference entre d'une part, une 
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distance entre le transducteur i du reseau d' imagerie et le 
point de focalisation qO, et d' autre part, une distance 
entre le transducteur i du reseau d' imagerie et le point de 
focalisation q. 

Dans le cas ou certains transducteurs d' indice v du 
reseau d' imagerie ne sont pas situes contre la couche 
aberratrice, il est en outre souhaitable de corriger les 
signaux de reference les signaux e' v (t,q) correspondants 
par repropagation numerique des transducteurs virtuels 
(situes contre la couche aberratrice 3) aux transducteurs 
reels (separes de ladite couche 3 par du gel 7 ou autre), 
de facon connue en soi, par le procede inverse de celui 
decrit ci-dessus en ce qui concerne les reponses 
impulsionnelles . 

15 0n fait ensuite emettre par les transducteurs du 

reseau d' imagerie, par iterations, les estimations obtenues 
des signaux de reference e'i(t,q), puis on capte avec les 
memes transducteurs les signaux, Si(t >q ), retrodif fuses par 
le milieu heterogene dissipatif. 

20 Puis ' on modifie pour .1' iteration suivante ces 

signaux de reference e'i(t,q) de la maniere suivante : 

ou les valeurs ti(q) et ai(q> sont un retard et un facteur 
d' amplitude correctifs, calculees pour maximiser un critere 
25 de coherence C entre lesdits signaux retrodif fuses . 

Avantageusement, le critere C pourra etre le 
suivant : 



< 

C = — 



OU 



• - Si{t>tl) = s i (t)®e l {t 1 q) , © representant 1/ operation de 

30 convolution, 

et .<> represente une moyenne temporelle. 
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Dans ce processus d' optimisation, les valeurs ^(q) 
peuvent etre calculees a chaque iteration de fagon a 
maximiser une fonction d' intercorrelafcion, pour des 
transducteurs voisins du reseau d'imagerie, des signaux 
5 gi(t,q) et gi+i(t,g) susmentionnes, et les valeurs cti(q) 
peuvent etre calculees de maniere a egaliser sur l'indice i 
1' amplitude maximale des fonctions g x (t,q). 

On optimise ainsi les signaux de reference 
e'i(t,q), i allant de 1 a n, de fagon qu'ils produisent un 
10 signal acoustique focalise de fagon precise sur le point de 
focalisation q situe dans le cerveau. Ce processus 
d' optimisation a deja ete explicite plus en detail par 
Mallart et al. (The Van Cittert-Zernike theorem in pulse 
echo measurements, J. Acoust. Soc. Am. 90(5), novembre 
15 1991, p. 2716-2727 ; Adaptative focusing in scattering 
media through sound speed inhomogeneities : the Van Cittert 
Zernike approach and focusing criterion, J. Acoust- Soc. 
Am.* 96(6), decembre 1994, p. 3721-3732). 

Lorsque cette optimisation est terminee pour un 
20 point de focalisation q, par exemple apres 2 ou 3 
iterations lorsque le critere C a atteint une valeur 
predetermines (notamment proche de 2/3), on passe au point 
de focalisation q+1 suivant, etc. 

Les signaux de reference e'i(t,q> ainsi obtenus 
25 sont memorises, par exemple dans les memoires Ml-Mn. 

Une fois l'etape d' apprentissage terminee, on peut 
notamment realiser des images echographiques du cerveau 2, 
eventuellement a un rythme eleve pouvant atteindre la 
cadence d'un echographe standard, par exemple 20 a 30" 
30 images par seconde. Pour realiser chacune de ces images, on 
procede comme suit, pour chaque point de focalisation q 
appartenant au champ a observer : 

on fait emettre respectivement par les 
transducteurs Ti du reseau d'imagerie, lesdits signaux de 
35 reference ei(t,q), 



WO 02/32316 PCT/FR01/03208 

23 

- puis on fait capter par lesdits transducteurs du 
reseau d'imagerie des signaux si(t) retrodif fuses par le 
milieu viscoelastique, 

- on convolue le signal de reference ei(t,q) emis 
5 par chaque transducteur du reseau d'imagerie avec le signal 

retrodif fuse si(t) capte par ce transducteur, 

puis on somme les produits de convolution ainsi 

obtenus . 

On notera que les differentes operations 
10 susmentionnees realisees au cours de l'etape 
d'apprentissage ou de l'etape d'imagerie peuvent etre soit 
programmers dans 1' unite centrale CPU, soit realisees en 
tout ou partie par des circuits specialises. 

Par ailleurs, on notera egalement que tous les 
15 transducteurs Ti, T r r pourraient etre utilises pour 
realiser les images echographiques du cerveau. Dans ce cas, 
le reseau d'imagerie serait confondu avec le reseau cible 
et. chacun de ces deux reseaux.. comprendrait tous -les 
transducteurs, le f onctionnement decrit ci-dessus 
20 s'appliquant alors mutatis mutandis. 
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RE VEND I CAT I ON S 

1. Procede non invasif de focalisation d'ondes 
acoustiques dans un milieu heterogene dissipatif (2,3) 
5 comprenant un milieu sensiblement homogene (2) entoure au 
moins partiellement par une couche aberratrice dissipative 
(3) qui genere des aberrations dans la propagation des 
ondes acoustiques, les ondes acoustiques etant emises 
depuis l'exterieur de la couche aberratrice (3) et 
10 focalisees dans le milieu sensiblement homogene (2), 

caracterise en ce qu'il comporte les etapes suivantes : 

(a) une etape initiale de positionnement au cours 
de laquelle on fixe un nombre t superieur a 2 de 
transducteurs acoustiques (Tl-Tn, T'l-T'm) dans des 

15 positions predetermines a l'exterieur de la couche 
aberratrice (3), ces transducteurs etant en contact avec 
ladite couche aberratrice et formant au moins : 

un reseau d' imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteurs, 

20 - et un reseau cible (T'l-T'm) qui regroupe un 

nombre m compris entre 1 et t desdits transducteurs, 

(b) une etape d' apprentissage comprenant elle-meme 
les sous etapes suivantes : 

(bl) une sous-etape d' apprentissage de 
25 focalisation du reseau d' imagerie sur le reseau cible, 
sous-etape au cours de laquelle : 

(bll) on determine des reponses 

impulsionnelles hri(t) du milieu heterogene dissipatif, 
respectivement entre chaque transducteur i du reseau 
30 d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice (3) en correspondance respective avec 
des transducteurs du reseau cible, ces reponses 
impulsionnelles etant memorisees sous forme numerique avec 
un certain echantillonnage temporel qui determine un nombre 
35 p de composantes frequent ielles de la reponse 
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impulsionnelle, de frequences respectives «k, i etant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d'imagerie, r etant un indice compris entre 1 et m 
qui designe un point de focalisation correspondant a un 
transducteur du reseau cible et k etant un indice compris 
entre 1 et p qui designe une composante f requentielle, 

(bl2) a partir de ces reponses 
impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un ensemble de n signaux de reference temporels 
e'i(t,r)., i variant entre 1 et n, tels que, si la paroi 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les differents transducteurs i du reseau d'imagerie 
15 genererait un signal predetermine focalise sur le point de 
focalisation r, 

(b2) une sous-etape de focalisation en un nombre 
R de- points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 
et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation lam correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(t,q) a faire emettre par les differents transducteurs i 
5 du reseau d'imagerie pour generer une impulsion focalisee 
sur ledit point de focalisation q, les signaux de reference 
e'i(t,q) etant determines pour chaque point de focalisation 
q en procedant coirane suit : 

(b21) une premiere estimation de e'i(t,q), 
pour i allant de 1 a n, est calculee a partir d'au moins un 
signal de reference e'i(t,q0), q0 etant 1 ' indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
ondes acoustiques dans le milieu sensiblement homogene (2), 
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(b22) on faifc emettre par les transducteurs 
ciu reseau d'imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e / i(t / q) / 
puis on capte avec les memes transducteurs des signaux 
5 Sj.(t,q) retrodif fuses par le milieu heterogene dissipatif, 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

Cf(0 a-M)>eM-rM)) 
ou les valeurs ai(q) et Ti(q) sont un facteur d' amplitude 

10 et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses, 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 

15 e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i(t,g), i 

20 etant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 

2. Procede selon la revendication 1, dans lequel au 
cours de la sous-etape (bll), lorsqu'au moins certains 
transducteurs (Tl-Tm, T'l-T'm) sont en contact avec un 

25 milieu homogene intermediaire lui-meme en contact avec la 
couche aberratrice, on corrige les reponses impulsionnelles 
hri(t) par repropagation numerique pour simuler des 
transducteurs situes directement au contact de la couche 
aberratrice. 

30 3. Procede selon la revendication 1 ou la 

revendication 2, dans lequel la sous-etape (bl2) comporte 
elle-meme les sous-etapes suivantes : 

(bl21) on determine p matrices de transfert 
H (cok) = [ Hri (cok) ] / i allant de 1 a n et r allant de 1 a m, 
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ou Hri(cok) est la valeur, a la frequence cok, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 
(bl22) on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
5 cible, n composantes Ei(cok,r), i variant entre 1 et n, 
telles que F (cok, r ) =H (cok) .E(cok,r), ou E(o)k,r> = [Ex (cok, r) ] 
est un vecteur a n composantes, F(cok,r) est un vecteur a m 
composantes Fl ((Dk/r), 1 variant entre 1 et m, ces m 
composan v tek (/ Fl (cok, r) correspondant a une focalisation 
10 "souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant a un transducteur du 
reseau cible, 

(bl23) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
15 cible, un vecteur de n signaux temporels e (t, r) = [ei (t , r) ) , 

/> 

i variant entre 1 et n, ou 0(/,r)=g£/(<afc t r).g i en notation 

complexe, ces signaux ei(t,r) etant adaptes pour que 
1' emission de ces derniers . respectivement par les 
diff§rents transducteurs i du reseau d'imagerie genere une 
10 impulsion acoustique focalisee sur le point de focalisation 
r du reseau cible, 

(bl24) une sous-etape de correction des 

aberrations generees par la couche aberratrice entre le 
milieu sensiblement homogene et chaque transducteur cible 
15 r, ces aberrations etant estimees sur la base des mesures 
precedemment effectuees, les aberrations ainsi estimees 
etant utilisees pour calculer lesdits signaux temporels de 
reference e'i(t,r) . 

4. Procede selon la revendication 3, dans lequel au 
0 cours de la sous-etape (b!22) on calcule p matrices 
H~ x (a>k), respectivement par regularisation et inversion des 
matrices de transfert H(cok>, et pour chaque transducteur r 
du reseau cible, on calcule le vecteur E(o)k,r> par la 
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formule : 

E(o)k,r)=H" i (o)k) .F{CDk,j) . 

5. Procede selon la revendication 3 ou la 
revendication 4, dans lequel au cours de I'etape (bl22), 

5 les composantes Fl(a)k,r) du vecteur F(o)k,r) correspondant 
a la repartition spatiale du champ desire a la frequence 
cok, sont egales a 0 pour l^r et egale a 1 pour 1-r. 

6. Procede selon l f une quelconque des 
revendications 3 a 5, dans lequel au cours de la sous-etape 

10 (bl24), on assimile la paroi aberratrice, au voisinage de 
chaque -point de focalisation r correspondant a un 
transducteur du reseau cible, a un filtre a reponse 
impulsionnelle finie, defini a chaque frequence cok par une 
amplitude Gr(cok) et une phase (^(cok), la sous-etape (bl24) 

15 comportant elle-meme les sous-etapes suivantes : 

(bl241) on calcule, pour chaque frequence 
cok, 1' amplitude Gr(cok) et la phase <t> r {CDk}, a partir soit 
des signaux ei(t,r), soit des vecteurs E(a>k, r) , 

(bl242) on calcule p matrices de transfert 

20 corrigees H' (o>k) = [H' ji(cok)], ou 

J Gj{CO k ) 

(b!243) on determine pour chaque 
transducteur r du reseau cible, n composantes E'i(cok,r), i 
variant entre 1 et n, telles que F (cok, r > =H r (ok) . E' (o>k, r ) , 

25 ou EMcak f r)= [E , i{o)k / r)] est un vecteur a n composantes, 
F(ook,r) est un vecteur a m composantes Fl(cok,r), 1 variant 
entre 1 et m, ces m composantes Fl(a)k,r) correspondant a 
une focalisation souhaitee des ondes acoustiques a la 
frequence cok sur le point de focalisation r correspondant 

30 a un transducteur du reseau cible, 

(bl244) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
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cible, un vecteur de n signaux temporels de reference 
e' (t,r) = [e'i(t,r) ), i variant entre 1 et n, ou 

e*t(t 1 r)=^Fi(a)k,r).g t ' en notation complexe. 

7. Procede selon la revendication 6, dans lequel au 
cours de la sous-etape (bl241>, on calcule 1' amplitude 
Gr(cok) et la phase <|> r (cok) comme suit : 



/=] 



^E(cQkj').Ei\cok,r) 



10 ou : 

Ei* est la valeur complexe conjuguee de Ei, 
et Ax(r0,r,i)=(d(r0,i)-d(r,i) )/c, d(r,i) etant la 
distance entre" le * transducteur i et ie point de 
focalisation r, et d(r0,i) etant la distance entre le 
15 transducteur i et un point de focalisation particulier rO. 

8. Procede selon la revendication 1 ou la 
revendication 2, dans lequel la sous-etape (b!2) comporte 
elle-meme les sous-etapes suivantes : 

(bl21) on determine p matrices de transfert 
20 H(a>k)-[Hri(©k) ], i allant de 1 a n et r allant de 1 a m, 
ou Hri(cok) est la valeur, a la frequence cok, de la 
transformer de Fourier de la reponse impulsionnelle hri(t), 
(bl22' ) on corrige les matrices de transfert 
H(ok) pour s'affranchir des aberrations generees par la 
25 paroi aberratrice au voisinage de ' chaque point de 
focalisation r, cette correction etant effectuee a partir 
des reponses impulsionnelles hri(t) precedemment 
determinees, et on obtient ainsi des matrices de transfert 
corrigees H' (ok) , 
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(bl23') on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, n composantes E'i(G)k,r), i variant entre 1 et n, 
telles que F(o>k,r)=H' (cok) . E' ((uk,r) , ou E' (cok,r) = 
5 [E'i(a>k,r)J est un vecteur a n composantes, F(cok, r) est un 
vecteur a m composantes Fl(a>k,r), 1 variant entre 1 et m, 
ces m composantes Fl(a>k,r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant a un 
10 transducteur du reseau cible, 

(bl24') on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels 
e ' (t,r) = [e'i(t,r) ] , i variant entre 1 et n, ou 
p 

15 e\{i,r)=^Ei{&ks). e iwkt en notation complexe, les signaux 
k=\ 

e'i(t,r) etant lesdits signaux de reference. 

9. Procede selon la revendication 8, dans lequel au 
cours de la sous-etape (bl23' } on calcule p matrices H'~ 
Mcok), respectivement par regularisation et inversion des 

20 matrices de transfert H' (cok) , et pour chaque transducteur 
r du reseau cible, on calcule le vecteur E' (cok, r) par la 
formule : 

E' (Q)k / r)=H'" :i (cDk) .F(cok, j) . 

10. Procede selon la revendication 8 ou la 
25 revendication 9, dans lequel au cours de 1 ' etape (bl23 7 ) , 

les composantes Fl(cok,r) du vecteur F(a)k,r} correspondant 
a la repartition spatiale du champ desire a la frequence 
cok, sont egales a 0 pour l*r et egale a 1 pour l=r. 

11. Procede selon 1 ' une quelconque des 
30 revendications 8 a 10, dans lequel au cours de la sous- 
etape (bl22')/ on assimile la parci aberratrice, au 
voisinage de chaque point de focalisation r correspondant a 
un transducteur du reseau cible, a un filtre a reponse 
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impulsionnelle finie, defini a chaque frequence cok par une 
amplitude Gr(cok) et une phase (|) r (a)k), la sous-etape (bl22' ) 
comportant elle-meme les sous-etapes suivanhes : 

(bl22'l) on calcule, pour chaque frequence 
cok, 1' amplitude Gr(cok) et la phase (^(rok), a partir des 
reponses impulsionnelles precedemment determinees, 

(bl22 f 2) on calcule p matrices de transfert 
corrigees H' (cok) « [H' ji (©k) ] , oil 

12. Procede selon la revendication 11, dans lequel 
au cours de la sous-etape (bl22'l)on calcule, pour chaque 
frequence cok, 1' amplitude Gr(cok) et la phase (|> r (cok), de la 
maniere suivante : 



1 5 Gr(0)k) 



_ i 



i n 

n^H^i{cok),H^j \cok) 



M™>)^ , ou : 

H*ri designe la valeur complexe conjuguee de Hri, 
et AT(rO,r,i)-(d(rO,i)-d(r,i))/c, d(r,i) etant la 
20 distance entre le transducteur i et le point de 
focalisation r, et d(r0,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO . 

13. Procede selon l'une quelconque des 
revendications precedentes, dans lequel au cours de 1'etape 
25 (c)., on fait suivre la sous-etape (cl) par les sous-etapes 
suivantes : 

(c2) on fait capter par lesdits transducteurs du 
reseau. d'imagerie des signaux s x (t) retrodif fuses par le 
milieu heterogene dissipatif, 
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(c3) on convolue le signal de reference emis par 
chaque transducteur du reseau d'imagerie avec le signal 
retrodiffuse capte par ce transducteur, 

(c4) puis on somme les produits de convolution 
5 ainsi obtenus, 

1'etape (c) etant renouvelee pour une pluralite de points 
situes dans le milieu sensiblement homogene . 

14. Procede selon l'une quelconque des 
revendications precedentes, dans lequel au cours de la 

10 sous-etape (b21) , la premiere estimation de chaque signal 
de reference est e'i(t,q) = e' i (ts+0i (q) , qO) pour chaque 
point de focalisation q, qO etant 1 ' indice d'un point de 
focalisation proche du point q pour lequel le signal de 
reference a deja ete determine, 9i(q> etant un retard egal 

15 a une valeur 5i(q)/c, ou c est la celerite moyenne des 
ondes acoustiques dans le milieu, et 5i (q) est egal a une 
difference entre d' une part, une distance entre le 
transducteur i du reseau d' imagerie et le point de 
focalisation qO, et d' autre part, une distance entre le 

20 transducteur i du reseau d'imagerie et le point de 
focalisation q, 

15. Procede selon l'une quelconque des 
revendications precedentes, dans lequel au cours de la 
sous-etape (b2), lorsqu'au moins certains transducteurs 

25 d' indice v du reseau d'imagerie ne sont pas directement au 
contact de la couche aberratrice, on corrige les signaux 
e'v(t,q) correspondants par repropagation numerique pour 
simuler des transducteurs places en contact direct avec la 
couche aberratrice . 

30 16. Procede selon l'une quelconque des 

revendications precedentes, dans lequel au cours de la 
sous-etape (b22), on recherche les valeurs a x (q) et x x (q) 
pour maximiser le critere de coherence C suivant : 
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/, -S < k*i( , - r /^)| 2 > 

£/('»(r) = */(0® ® representant 1' operation de 

convolution, 

et <> represente une rnoyenne temporelle. 
5 17. Procede selon la revendication 16, dans lequel 

au cours de la sous-etape (b22>, les valeurs ti(q) sont 
calculees en maximisant une fonction d' intercorrelation, 
pour des transducteurs voisins du reseau d'imagerie, des 
signaux gi(t,q) et gi+i ( t , q ) . 
10 18 • Procede selon la revendication 16 ou la 

revendication 17, dans lequel au cours de la sous-etape 
(b22), les valeurs <Xi(q) sont calculees de maniere a 
egaliser sur l'indice i 1' amplitude maximale des fonctions 
gi (t/q) . 

I 5 19 • Procede selon l f une quelconque des 

revehdications precedentes, dans lequel la. sous-etape (b22> 
relative a chaque point de focalisation q est realisee 
iminediatement apres la sous-etape (b21> relative au meme 
point de focalisation q. 

20 20. Procede selon l'une quelconque des 

revendications precedentes, dans lequel le milieu 
heterogene dissipatif est constitue par le cerveau entoure 
par le crane. 

21. Procede selon l'une quelconque des 
25 revendications precedentes, dans lequel : 

- soit le reseau d'imagerie et le reseau cible 
sont deux reseaux distincts disposes de part et d 1 autre du 
milieu heterogene dissipatif, 

- soit tous les transducteurs appartiennent a la 
0 fois au reseau d'imagerie et au reseau cible. 

22. Procede selon l'une quelconque des 
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revendications precedentes, dans lequel les ondes 
acoustiques sont des ondes ultrasonores . 

23. Dispositif (1) adapte pour mettre en ceuvre un 
procede selon l'une quelconque des revendications 
5 precedentes, ce dispositif comportant un nombre t superieur 
a 2 de transducteurs acoustiques (Tl-Tn, T'l-T'm) destines 
a etre fixes dans des positions predetermines a 
I'exterieur de la couche aberratrice (3), ces transducteurs 
etant commandes par au moins une unite centrale 

10 electronique (CPU) et forraant au moins : 

un reseau d'imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteurs, 

- et un reseau cible (T'l-T'm) qui regroupe un 
nombre m compris entre 1 et t desdits transducteurs, 

15 1' unite centrale electronique etant adaptee pour suivre les 
etapes suivantes : 

(b) une etape d' apprentissage comprenant elle-meme 
les sous etapes suivantes : 

(bl) une sous-etape d' apprentissage de 

20 focalisation du reseau d'imagerie sur le reseau cible, 
etape au cours de laquelle : 

(bll) on determine des reponses 

impulsionnelles hri(t) du milieu heterogene dissipatif, 
respectivement entre chaque transducteur i du reseau 

25 d'imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible, ces reponses impulsionnelles 
etant memorisees sous forme numerique avec un certain 
echantillonnage temporel qui determine un nombre p de 

30 composantes f requentielles de la reponse impulsionnelle, de 
frequences respectives o)k, i etant un indice compris entre 
1 et n qui designe un transducteur du reseau d'imagerie, r 
etant un indice compris entre 1 et m qui designe un point 
de focalisation correspondant a un transducteur du reseau 

35 cible et k etant un indice compris entre 1 et p qui designe 
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une composante frequentielle, 

(bl2) a partir de ces reponses 
impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la paroi 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les differents transducteurs i du reseau d'imagerie 
genererait une impulsion acoustique focalisee sur le point 
de focalisation r, 

(b2) une sous-etape de focalisation en un nombre 
R de points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 
15 et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation lam correspondant aux 
transducteurs du reseau oible, des signaux de reference 
e'i(t,q) a faire emettre par les differents transducteurs i 
20 du reseau d'imagerie pour generer un signal predetermine 
focalise sur ledit point de focalisation q, les signaux de 
reference e'i(t,q) etant determines pour chaque point de 
focalisation q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 
25 pour i allant de 1 a 4, est calculee a partir d'au moins un 
signal de reference e'i(t,q0), q0 etant 1 ' indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
30 ondes acoustiques dans le milieu sensiblement homogene (2), 

(b22) on fait emettre par les transducteurs 
du reseau d'imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q), 
puis on capte avec les memes transducteurs des signaux, 
35 s A (t,q) retrodif fuses par le milieu heterogene dissipatif, 
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puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

e](t) -> a l (c I ).e i {t-T i ((/)) 
ou les valeurs cti(q) et ii(q) sont un facteur d' amplitude 
5 et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses , 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 
10 e'i(t,q), au moins pour q compris entre m+1 et m+R, 

'(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d ! imagerie, lesdits signaux de reference e'i(t,q), i 
15 etant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 
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